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Abstract 

A series of simulation studies (Ono & Ikegami 1999; 
2001) show that a proto cell spontaneously emerges from 
a chemical soup by acquiring membrane structures. In 
2-dimensional space, the emergence of proto cells is fol- 
lowed by the reproduction of cells. A major unsolved 
problem is the evolution of proto cells; how the proto 
cells evolve into modern cells with higher functionali- 
ties. Here we examine, as the first step, the evolution 
of catalysts within the proto-cells. Catalytic chemicals 
have different catalytic activity in generating membrane 
chemicals. We show that cells with higher activity of 
membrane production evolve through cellular selection. 



Introduction 

It is widely accepted that the origin of life was a set 
of molecules that catalyzed the reproduction of each 
other. However, when we consider the evolution of 
such primitive chemical systems, the compartment aliza- 
tion of molecules is indispensable for establishing the 
co-evolution of cooperative catalytic reactions and pro- 
tecting them from parasites that would spoil the evolu- 
tion (Szathmary & Maynard Smith 1997). Though it is 
difficult to know about the structure of primitive cells 
because there remain few physical records of the earli- 
est living cells, there have been various theoretical ap- 
proaches to understanding the emergence and evolution 
of proto-cell systems. 

Ganti proposed a minimal model of primitive self- 
maintaining cells named "chemoton" (Ganti 1975; 1997). 
It is composed of (1) a metabolic system of autocatalytic 
molecules, (2) self-replicating molecules that inherit ge- 
netic information and (3) a self-organizing membrane 
molecule to enclose the whole system. This system main- 
tains itself by consuming resources and discharging waste 
into the environment. It can be easily imagined that if 
the reproduction of the cell is appended to this system, 
it would be a primitive unit of Darwinian selection and 
evolve into more stable structures. 

It should be stressed that a cell defines itself as an 
individual by producing a membrane that distinguishes 
itself from the outside. Mat ur ana and Varela pointed 



out that this is a unique feature of living organisms, 
and named it "autopoiesis" (Maturana & Varela 1980). 
To demonstrate the self-maintenance of an autopoietic 
structure, abstract computational models of an autopoi- 
etic cell based on a Cellular Automaton were proposed 
(originally by Varela (Varela, Maturana, & Uribe 1974), 
and re-implemented by Zeleny (Zeleny 1977) and by Mc- 
Mullin (McMullin & Varela 1997)). Breyer and Mc- 
Caskill introduced the metabolism of a catalyst into 
this model (Breyer, Ackermann, & McCaskill 1998). 
It was also shown that an autopoietic proto-cell can 
reproduce itself automatically (Ono & Ikegami 1999; 
2001). Speroni di Fenizio and Dittrich proposed another 
approach to represent proto-cells that are embedded in 
a triangular planar graph (Speroni di Fenizio, Dittrich, 
& Banzhaf 2001). 

The remaining question is "How was the first cell orga- 
nized?" Answering this question will give us the first step 
in understanding the emergence of higher order struc- 
tures in life's evolution. This paper consists of two parts. 
The computational algorithm of the model is explained 
in detail in the first part. We introduce a Lattice Artifi- 
cial Chemistry (LAC) model that simulates the chemical 
reactions and spatial interactions of abstract chemicals. 
In the second part, an emergence of a proto-cell from a 
non-organized initial state, its reproduction and the se- 
lection of inside catalysts through the cell reproduction 
are reported in order. 

Lattice Artificial Chemistry 

Our model is based on discrete and stochastic dynamics, 
which is extended from a lattice-gas model. Chemicals 
are represented by particles on reaction sites that are 
arranged as a two-dimensional triangular lattice. Note 
that any number of particles can be placed on a single 
site. The vector n(x) = (ni(x), n2(x), . . . , nm(x)) gives 
the number of each type of particles on the site x. Ni 
gives the total amount of z-th particles in the system. 

Chemical reactions are expressed by the probabilistic 
transition of particle types. The diffusion of chemicals 
is expressed by random walks of particles on the sites. 
These transition probabilities are given as the products 



of the associated rate coefficients and the fohowing func- 
tion of the potential change AE^ 

/(^^) = ^p^y (1) 

where /3 represents the inverse of the product of 
the Boltzmann constant and temperature (note that, 
f{AE)/f{-AE) = e-^^^). In the simulations reported 
hereafter, the value of P is normalized and fixed to 1. 

Hydrophobic Interaction 

The probabilities of random walks of particles are biased 
according to the gradient of the potential ^(x) which is 
given by summing up the interaction from all particles 
in the same and adjacent sites. The probability p with 
which a particle i moves from a site x to x' is calculated 
as follows, 

|x'-x|<l j 

Pi(x^x') = Dz/, /(vl/,(x') - \E',(x)) (3) 

where ^i(x) denotes the potential of particle i in the 
site ^^Dif^ denotes the diffusion coefficient of particle i, 
and ipij{dx) denotes the interaction on particle i from 
particle j. Diffusion coefficients depend on the species 
of the particles. Autocatalytic and membrane particles 
are assumed to be larger molecules so that their diffu- 
sion coefficients are smaller than those of other particles 
(Di/^. = Difu = 0.003, Dif^,^^^,, = 0.01). 

To simulate the formation of membranes, we intro- 
duce hydrophobic interactions between particles. First, 
the particles are grouped into three classes: hydrophilic, 
hydrophobic and neutral. In general, all particles repel 
each other, but repulsion between hydrophilic and hy- 
drophobic particles is much stronger than that between 
other particles so that phase separation between different 
classes of particles takes place. On the other hand, neu- 
tral particles do not repel other particles very strongly 
so that they can diffuse more freely. 

Next, we assume that hydrophobic particles M are 
anisotropic. Namely, the repulsion around a particle M 
depends on its orientation and the configuration of the 
particles as illustrated in Fig. 1. There are specific direc- 
tions in which the repulsion becomes strong, while the 
repulsion becomes weak in the other directions Taking 
its symmetry into account, a particle M can rotate to 
six different orientations (M^, M±^/^, M±^/^, M^/^) 
stochastically according to the gradient of the potential 
as follows, 

V^,^y{-s) = Rotf{^^y{^)-^^k{^)) (4) 

where Rot = 0.01 denotes the rotation coefficient. The 
repulsion between two particles M becomes strong when 
their orientations are different, so that they tend to align 
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Figure 1: Illustration of repulsion around a particle M. 
The honeycomb cells represent the lattice sites. The 
depth of gray shade corresponds to the magnitude of 
repulsion against a hydrophilic particle on the site from 
the particle M on the center site. The repulsion becomes 
stronger on the dark gray sites than on the light gray 
sites. 

in the same orientation. According to these interactions, 
particles M gather together to form stretched clusters. 
We call these stretched structures of particle M "mem- 
brane". Though the characters of membranes, such as 
flexibility, depends on these values, the formation of 
membranes can be observed in a wide range of param- 
eters. The detail values of repulsion ipij^d^) which are 
arbitrary chosen for the following experiments are listed 
in Tables la and lb. 



Table 1: a. Repulsion between isotropic particles. 







position 


particles 


= 


dr = 1 


hydrophilic 


hydrophilic 


0.0100 


0.0033 




neutral 


0.0010 


0.0003 


neutral 


neutral 


0.0010 


0.0003 



Chemical Reaction 

We introduce a simple metabolic system of autocatalytic 
particles. Consider that there are various species of self- 
replicating particles, and some of them have the ability 
to catalyze the production of membrane particles. Re- 
sources of these particles are supplied from some external 
source homogeneously. 

Figure 2 illustrates the reaction paths. The probabili- 
ties of chemical reactions depend on the enthalpy change 
along with the transition, as follows, 

AHij = AHj - AHi (5) 
p,^,(x) = fc,,(x)/(Ai7,,+^,(x)-^,(x)) (6) 

where AHij denotes the enthalpy change that is given 
by the difference in the formation enthalpy listed in Ta- 
ble 2. 



Table l.b Repulsion between hydrophobic and other particles 



particles 


position 


dr = 


dr = 1 


|9 = 0,7r 


= 7r/3, -27r/3 


6 = 27r/3, -7r/3 


MO 


hydrophilic 


0.2000 


0.1600 


0.0200 


0.0200 


neutral 


0.0010 


0.0008 


0.0001 


0.0001 




0.0100 


0.0033 


0.0033 


0.0033 




0.0777 


0.0259 


0.0259 


0.0259 




0.1433 


0.0477 


0.0477 


0.0477 




0.2100 


0.0700 


0.0700 


0.0700 




hydrophilic 


0.2000 


0.1000 


0.1000 


0.0000 


neutral 


0.0010 


0.0005 


0.0005 


0.0000 




0.0100 


0.0033 


0.0033 


0.0033 




0.2100 


0.0700 


0.0700 


0.0700 



kij(yi) denotes the coefficient of reaction i ^ j that 
may depend on the number of catalysts on the site. Note 
that the effects of the interactive potential, namely, the 
effects of hydrophilic/hydrophobic environments are also 
taken into account here, therefore, for example, it be- 
comes more difficult to synthesize a hydrophilic particle 
in a hydrophobic environment. 
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Figure 2: Schematic drawings of reaction paths. An 
autocatalyst (Ai) catalyzes the reproduction of another 
particle Aj from a resource particle (Xa) that has a 
higher chemical energy. It also catalyzes the produc- 
tion of a membrane particle (M) from another resource 
particle (Xm)- All particles decay into waste particles 
(Ya and Ym, respectively), spontaneously, however, an 
external energy supply recycles Ya and Ym into Xa 
and Xm, respectively. The number of total particles is 
preserved. 



tion of another particle Aj using it as a template and 
consuming a resource particle (Xa). 



Ai 
AiAj^ 



Xa 



AiAj 
AiAj 



(7) 
(8) 



There is a probability of mutation /i with which a particle 
Ai mutates to Ai±i when it is reproduced. Assuming 
that the rate of the ffist reaction is much faster than that 
of the second one, the rate coefficients between Xa and 
Aj can be given as follows. 



. (x) = /inA,._, (x) + (1 - 2/i)nA,. (x) + /in^^.^i (x) (9) 



:(X) 



/cA + CAn:4^(x)^n^,(x) (10) 



where (x) denotes the number of particles Ai on the 
site X, and kA denotes the rate of spontaneous reaction. 
Note that all autocatalysts share a common catalytic 
activity Ca and catalyze the replication of each other 
equally. 

An autocatalytic particle also catalyzes the produc- 
tion of a membrane particle (M) from a resource (Xm)- 
The activity {Cmi) depends on the species. The cat- 
alytic activity of each species Ai is given by the follow- 
ing equation, namely, the activity of particle Ai is i-th 
times larger than that of particle Ai, so that the rate 
coefficients between Xm and M are given as follows. 

Cm xi (11) 
kxM^Mi^) A;m + ^C'M,nA,(x), (12) 



Table 2: Formation enthalpy 



particle 


Xa, Xm 


Ai 


M 


Ya, Ym 


AH, 


12.0 


6.0 


4.0 


0.0 



There are ten species of autocatalytic particles ( Aq . . . 
Ag). An autocatalytic particle Ai catalyzes the replica- 



where Cm is a given constant, and kM denotes the rate 
of spontaneous reaction. 

These particles naturally decay into waste particles 
(Ya and Ym, respectively) at a constant rate ky- How- 
ever, we introduce an external source that supplies re- 
sources. To preserve the total number of particles, the 
resource supply is expressed by the exchange from waste 



to resource particles. Thus the transition coefficients are 
given as fohows, 



kxA^YA ^Xm^Ym 
^Ya^Xa ^Ym^Xm 



ky^Sx- 



(13) 
(14) 
(15) 



Due to the term Sx^ the whole system is kept in a non- 
equilibrium state. The last particle (W) represents wa- 
ter that does not change into other particles. We assume 
that water and autocatalytic particles are hydrophilic 
particles that are repelled by membranes, and resource 
and waste particles are neutral particles which can dif- 
fuse through membranes. 

The rate coefficients of spontaneous reactions are kA = 
ku = 1.0 X 10"^, ky = 1.0 X 10"^. The coefficients 
of catalytic activity are Ca = 1.0 x 10~^ and Cm = 
1.0 X 10"^. The mutation rate is /i = 1.0 x 10"^^. The 
rate of resource supply is given a constant Sx = 16. 

The simulation is based on a Metropolice method. At 
each iteration, the following steps are repeated, 

1. Calculate the potential of each particle. 

2. Calculate the probabilities of diffusion, rotation and 
chemical transition according to the potential differ- 
ence. 

3. Change the state of particles according to the proba- 
bilities synchronously. 

In the initial state, the particles are placed randomly. 
There are 30 particles on a site on average, and the mean 
numbers of particles on a site are listed in Table 3. There 
is a sufficient number of resource particles and supplies 
to sustain metabolism. The average production rate of 
membranes is set very low at ffist. Catalysts with higher 
activity only emerge through random mutations. The 
reaction sites are arranged as a 64 x 64 triangular lattice 
whose boundaries are periodic. 



Table 3: Mean numbers on a site in the initial state 



particle 


Ao 


Ai 


A2 


A3 


A4. 


.Ag 


Hi 


1.6 


1.2 


0.8 


0.4 


0.0 


particle 


Xa 


Ya 


M 


Xm 


Ym 


w 


Hi 


2.0 


2.0 


0.0 


5.0 


5.0 


10.0 



Simulation Results 

The evolution of this system is roughly divided into three 
characteristic stages: (1) Chemical evolution, (2) Emer- 
gence of proto-cells, and (3) Cellular evolution. 



Chemical Evolution 

Fig. 3(1) shows the initial configuration. Before cellular 
selection starts, the chemical evolution simply depends 
on the reproduction rate of each species. In this model, 
because they share the same reproduction rate, the evo- 
lution is mostly driven by mutations and random fluctu- 
ation. At ffist, the largest part of the autocatalytic par- 
ticles is Ao which does not produce membrane particles. 
Table 4 shows a profile of the population after 30,000 
iterations for a single run. However, as the populations 
of other species increase, small pieces of membranes are 
gradually formed. (Fig. 3(2)). 





(1) Initial state 



(2) After 30,000 iterations 



Figure 3: Chemical evolution. The white regions are 
dominated by particle M. The depth of gray shade rep- 
resents the total population of the autocatalysts {Y^ Ai). 
The black regions are dominated by particle W. Re- 
source and waste particles are not displayed in the fig- 
ures. Pieces of membranes are produced by the catalysts 
which emerged through mutations. 



Table 4: A profile of the population of particles after 



particle 


Ao 


Ai 


A2 


A3 


A4 


rii 


1.00 


0.60 


0.42 


0.20 


0.05 


particle 


A5 


Ae 


At 


As 


Ag 


Hi 


0.07 


0.08 


0.05 


0.04 


0.06 


particle 


Xa 


Ya 


M 


Xm 


Ym 


Hi 


3.00 


4.43 


1.57 


3.68 


4.77 



Emergence of Proto-cells 

Once membranes are formed, they begin to restrict the 
diffusion of catalysts. Thus, membranes can keep the 
local population and also their reaction rate high. As 
resource particles are consumed faster in such regions, 
resource particles diffuse into these regions according to 
the gradient of the population. It increases their reaction 
rate more. Due to this osmotic competition for resources, 
a small difference in the population of autocatalysts be- 
tween the two sides of the membrane becomes larger. 




(1) 60,000 iterations (2) 120,000 iterations 




(3) 180,000 iterations (4) 240,000 iterations 



Figure 4: Emergence of Proto-cell structures. As the 
membranes grow, competition for resources between re- 
gions separated by membranes takes place. The regions 
differentiate into two states. In some regions that are en- 
closed by membranes, the density of autocatalysts stays 
high. In the other regions, their density becomes almost 
zero. 



When the density of resources becomes too low in 
some regions, autocatalysts are no longer able to sustain 
their replication. Autocatalysts and membrane particles 
in these regions decay. At last, most regions become in- 
active, namely, filled with only resource and waste parti- 
cles (and water), while there remain some active regions 
in which autocatalysts keep reproducing. These struc- 
tures maintain themselves autopoietically, namely, the 
autocatalysts inside them reproduce themselves and me- 
tabolize the membrane particles to maintain their mem- 
branes. We call this structure a "proto-cell" hereafter. 

Cellular Reproduction 

A proto-cell can not only maintain itself but can repro- 
duce itself. Figure 5 shows snapshots of the reproduc- 
tion process. A proto-cell structure grows in size by as- 
similating resource particles from neighboring regions. 
As it grows, it comes to produce more membrane parti- 
cles than it needs to maintain its membrane. When it 
reaches a certain size, surplus membrane particles begin 
to form another membrane inside the cell. This divides 
the mother cell into a few daughter cells. The daughter 
cells can continue to grow and reproduce recursively. 

On the other hand, if a cell fails to maintain its mem- 
brane, the density of catalysts lowers quickly due to dif- 
fusion through the defect of membrane. When the den- 




(1) 300,000 iterations (2) 309,000 iterations 




(3) 318,000 iterations (4) 327,000 iterations 



Figure 5: Reproduction of proto-cell. Snapshots from 
300,000 iterations to 327,000 iterations. The proto-cell 
indicated by the white arrow grows gradually in size. 
When it becomes too large, another membrane appears 
inside it. At last, it divides the cell into daughter cells. 
On the other hand, the cell indicated by the gray arrows 
could not produce enough membrane particles to keep 
its membrane. Catalysts in the cell diffuse away through 
the defect of the membrane. 



sity becomes too low, the catalysts can not sustain the 
metabolism any longer, and the whole structure finally 
disappears. This mutual dependency between membrane 
and metabolic system is essential to the evolution of 
proto-cells. 

Selection of Catalysts Through Cellular 
Reproduction 

Figure 6 shows the evolution of the population of cat- 
alysts. If there is no spatial structure, the populations 
of catalysts just randomly drift around the equilibrium 
where all populations are equal. However, once the 
differentiation into active and inactive regions has ad- 
vanced, it becomes difficult to sustain metabolism with- 
out membranes, because the density of autocatalysts 
lowers fast. 

Note that when a proto-cell divides itself, the composi- 
tion of the population of replicators contained in the cell 
is roughly inherited by the daughter cells. This implies 
that catalysts that can sustain a proto-cell more stably 
are selected regarding a proto-cell as a new unit of Dar- 
winian evolution. Therefore, due to cellular selection. 



the populations of catalysts are biased toward higher 
membrane production activity against the random drift. 
The populations of catalysts and their relative amount 
{pi = rTi/^rTj x 100) after 600,000 iterations for the 
same run as the previous one are listed in Table 5. 
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Figure 6: Evolution of the population of autocatalysts. 
At first, the population of Aq, which does not pro- 
duce membrane particles, is the largest. However, other 
species arise through mutation, and after the formation 
of proto-cells, species which can produce membranes 
have an advantage over Aq (see the thick lines which 
denote species A7 - Ag). 



Table 5: A profile of the populations of catalysts after 
600,000 iterations. 



particle 


Ao 


Al 


A2 


A3 


A4 




0.24 


0.36 


0.28 


0.30 


0.23 


P^{%) 


7.0 


10.5 


8.1 


8.7 


6.7 


particle 


A5 


Ae 


At 


As 


Ag 


Hi 


0.31 


0.31 


0.46 


0.43 


0.52 


P^{%) 


9.0 


9.0 


13.4 


12.5 


15.1 



To make the effects of the cellular selection more clear, 
we investigated evolution under a lower resource sup- 
ply. The initial configuration and the parameters are the 
same as those of the previous run. After the formation 
of the proto-cell structures (after 300,000 iterations), the 
rate of resource supply Sx was decreased to Sx x 0.75. 
Because the density of resource particles decreases, cat- 
alysts have to keep their density high to maintain a suf- 
ficient reproduction rate more tightly, and a cell that 
fails to maintain its membrane disappears faster. The 
pressure to acquire higher membrane production activ- 
ity becomes stronger. Table 6 shows the populations 
after 900,000 iterations for another run. The dominance 
of catalysts with higher activity is clearer. 

To see the detailed process of cellular selec- 



Table 6: A profile of populations after 900,000 iterations 
under the lower resource supply. 



particle 


Ao 


Al 


A2 


A3 


A4 


Hi 


0.19 


0.10 


0.11 


0.16 


0.13 


P^{%) 


8.6 


4.5 


5.0 


7.3 


5.9 


particle 


A5 


Ae 


At 


As 


Ag 


fii 


0.18 


0.25 


0.20 


0.42 


0.46 


P^{%) 


8.2 


11.4 


9.1 


19.1 


20.9 



tion, the evolution of the mean activity {Cmi = 
^ Cmi I ) of membrane production is shown 

in Fig. 7. Within a cell, because competition among the 
autocatalysts is neutral, the mean activity of each cell 
fluctuates randomly. But among the cellular assembly, 
a proto-cell in which membrane production activity is 
low becomes extinct more often. As a result, cells with 
higher catalytic activity outperform the lower ones and 
the total average of gradually increases. 



Total average 
Average witnin a cell • 
Cell inheritance 
Extinction X 




100 200 300 400 500 600 700 



Figure 7: Evolution of the mean activity of membrane 
production (CmJ- A cell with lower membrane produc- 
tion breaks more often because its membrane is weaker. 

Though there often appear cells with very high ac- 
tivity of membrane production, the Cmi of these cells 
soon drops again. These drops are caused by the inva- 
sion of "parasitic" catalysts. Namely, the evolved cells 
are dominated by catalysts that have higher activity and 
produce enough membrane particles, however, the cat- 
alysts with lower activity, e.g., Ao, can always emerge 
through mutation and increase through the random fluc- 
tuation within these cells. Snapshots of the process of 
cellular selection are shown in Fig. 8. At 780,000 itera- 
tions, the cells indicated by the white arrows are deeply 
infected by parasites. It becomes difficult for these cells 
to maintain their membranes and they disappear before 
900,000 iterations. The cells at 900,000 iterations are 
produced from survived cells that have higher catalytic 
activity, but, there are cells that are newly invaded by 



parasites (indicated by the gray arrow). This result in- 
dicates that these proto-cehs have hmited "hfe spans". 
A proto-ceh has to keep dividing itself to escape from 
parasites, otherwise, the parasitic catalysts increase in it 
before long. 




(1) 780,000 iterations (2) 900,000 iterations 

Figure 8: Invasion of parasitic catalysts. The white sites 
are where the mean catalytic activity is low. 



As the cellular evolution proceeds, it is observed that 
the size of the cells becomes smaller. Figure 9 shows 
a histogram of cell sizes (namely, the number of sites 
within the membranes) for the same run, after 300,000, 
600,000 and 900,000 iterations. The mean sizes of cells 
at each iteration are 103.8, 82.7, and 66.2, respectively. 
When the mean activity of the membrane production of 
a cell increases, it can maintain its membrane with fewer 
catalysts, and it can divide itself faster. However, it is 
sometimes observed that a cell whose size is too small 
fails to maintain itself. When a cell produces too many 
membrane particles or the size of cell become too small, 
the rate of replication of the autocatalysts decreases, be- 
cause the hydrophobic environment suppress the synthe- 
sis of hydrophilic particles. A certain optimum size of 
the cell will be achieved through the selection of cata- 
lysts. This expectation is supported by the result that, 
in the evolved state, the dispersion of cell sizes becomes 
smaller. 

Conclusion and Discussion 

In this article, we have presented a model for the evo- 
lution from molecular to cellular reproduction. Starting 
from a homogeneous, random initial state, we demon- 
strated the emergence of proto-cells. This goes through 
three stages: (1) Metabolic cycles that produce mem- 
brane molecules can arise through pre-cellular chemical 
evolution. (2) Proto-cell structures, i.e., self-maintaining 
structures that maintain their own membranes by them- 
selves emerge spontaneously. (3) A proto-cell structure 
reproduces itself. Because the molecules contained inside 
a proto-cell can pass genetic information into daughter 
cells, it can be regarded as a unit of Darwinian evolu- 
tion. Cells that can maintain themselves more stably 
evolve through cellular selection. 
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Figure 9: The distribution of cell sizes after 300,000, 
600,000 and 900,000 iterations. The standard deviations 
of the sizes are 47.52, 26. 58, and 24.26, respectively. 



Excess production of membrane particles generates a 
hydrophobic environment, which is less optimal for the 
self-reproduction of the particles. Nevertheless, particles 
with higher catalytic ability replicate more via cellular 
level selection. As a result, the average value of catalyza- 
tion with membrane formation becomes higher than that 
without membrane formation. We insist here that the 
emergence of cellular structures produces a new rugged 
evolutionary landscape on the particle level. 

The entire behavior is insensitive to parameter val- 
ues whenever we have the conditions of, 1) a membrane 
formation and 2) transportation of resource particles 
through membranes. This robust behavior gives an ad- 
vantage to our model over other cellular formation sys- 
tems such as Grey-Scott, whose spatio-temporal pattern 
is sensitive to reaction and diffusion rates. In our model, 
the size of cells and replication rates change but their 
qualitative behaviors never change. Instead, our model 
behavior is sensitive to the form of repulsion potential 
of the membrane particles. In this study, we chose a set 
of values which provide flexible membranes to make the 
organization of proto-cells easier. We take the advan- 
tage of this sensitivity to the potential form in order to 
study the evolutionary dynamics of these characteristics 
of membranes of proto-cells. 

In the simulations reported here, the resources supply 
was homogeneous. And the pressure to evolve membrane 
production became clear when the resource supply rate 
was lowered. These results suggest that, if the supply of 
resource is not homogeneous, namely, if there is a small 
region where the resource supply is plenty enough to 
sustain metabolism, but the supply is rather poor in the 
other regions, pre-cellular evolution will take place where 
the resource supply is high enough, and once the proto- 
cell structures are acquired, they can migrate where the 
resource supply is lower, which will promote the further 



evolution of cellular structures. 

Though the model introduced in this article is sim- 
ple and abstract, we are now extending our model to 
implement more complex metabolic reactions that can 
produce more diverse membrane particles and various 
cell types. Our next objective is to investigate the co- 
evolution of metabolic systems and cellular structures. 
Because in our model, the birth and death of cells are 
all actualized through the elemental interactions of par- 
ticles, without any ad- hoc rules, the way in which the 
proto-cell is organized itself can evolve. In this sense, it 
can be regarded as a basic model for the organization 
of a dynamical hierarchy. We expect that higher order 
structures (e.g., cell differentiation and cooperative in- 
teraction between cells) are yielded. 

Furthermore, the dynamics of our model is based on 
a local equilibrium system. A quantitative analysis of 
this model may bring insight to the evolution of prim- 
itive cells as a non-equilibrium system from a thermo- 
dynamical aspect. Along these lines, computational 
models of artificial chemistry will provide useful tools 
for understanding the earliest evolution of life. 
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